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Nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) of Cu have been measured in a cou- 
pled spin dimers and chains compound Cu2Fe2Ge40i3. Cu NQR has also been measured in an isostructural material 
Cu2Sc2Ge40i3 including only spin dimers. Comparison of the temperature dependence of the *^Cu nuclear spin-lattice 
relaxation rate between the two compounds reveals that the Fe chains in Cu2Fe2Ge40i3 do not change a spin gap energy 
of the Cu dimers from that in Cu2Sc2Ge40i3, contributing additionally to the relaxation rate at the Cu site. A modestly 
large internal field of 3.39 T was observed at the Cu site in the antiferromagnetic state of Cu2Fe2Ge40i3 at 4.2 K, which 
is partly because of quantum reduction of the ordered moment of a Cu atom. The internal field and the ordered moment 
of Cu are noncoUinear due to large anisotropy of the hyperfine interaction at the Cu site. A model analysis of the internal 
field based on the fourfold planar coordination of Cu suggests that a 3rf hole of the Cu"^ ion is mainly in the d(x^ - y^) 
orbital state. 

KEYWORDS: NMR, NQR, Cu2Fe2Ge40i3, Cu2Sc2Ge40i3, spin dimer, spin chain, nuclear spin-lattice relaxation 
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1. Introduction 

Effects of quantum fluctuation on the magnetic properties 
of low-dimensional spin systems appear in a variety of ways 
depending on dimensionality of a system, a spin quantum 
number, geometry of interactions, and so on. Systems like 
spin dimers, ladders, integer spin chains have a disordered 
ground state with a gap in the spin excitation spectrum, be- 
ing characterized by a short range spin correlation. This state 
is robust for external perturbations such as three-dimensional 
interactions and an addition of spin defects, and is not easily 
transformed into a state with long-range magnetic order In 
contrast, systems with gapless excitations such as half-integer 
spin chains are more easily driven to have long-range order, 
reflecting a quantum critical nature of the ground state which 
has a divergently large spin correlation length. 

Recently, fascinating model systems combining both types 
of spin networks have been realized experimentally and have 
attracted much attention expecting exotic phases and transi- 
tions between them.'^ Cu2Fe2Ge40i3 is one of the members 
of such a new family of materials, comprising of quantum 
spin dimers of Cu^+ ions {S - 1/2) and classical spin chains of 
Fe-^"^ ions {S - 5/2). The crystal structure of Cu2Fe2Ge40i3 is 
shown in Fig. 1. Cu2Fe2Ge40i3 crystallizes in the monoclinic 
space group P2\lm with the unit cell dimensions a - 12.101 
A, /7 = 8.497 A, c = 4.869 A, and yS = 96.131° at room tem- 
perature.^'^ The structure consists of three groups of metal- 
oxygen polyhedra: crankshaft-shaped chains of edge-sharing 
FeOe octahedra running in the b direction, dimers of edge- 
sharing Cu04 squares bridging the FeOe chains along the a 
axis, and Ge04 tetrahedra separating the chain-dimer blocks 
stacked in the c direction. It can also be viewed as being build 
from blocks of the famous spin-Peierls compound CuGeOs^ 
and the Fe oxide chains. 

Cu2Fe2Ge40i3 orders antiferromagnetically below the 
Neel temperature - 39 K.^ The magnetic structure is 



roughly collinear with slight canting of ordered moments be- 
tween the Cu and Fe sublattices. The ordered moments of 
Cu and Fe are estimated to be 0.38 and 3.62 jUb at 1.5 K, 
respectively. There are large and moderate reductions of the 
Cu and Fe moments, indicating strong quantum fluctuations. 
In the paramagnetic state, a temperature dependence of the 
bulk magnetic susceptibility is roughly reproduced by a sum 
of contributions of isolated S - 5/2 chains and S - 1/2 dimers. 
This has rendered support to a naive picture of Cu2Fe2Ge40i3 
being a bicomponent magnet with weakly coupled spin chains 
and dimers. 

The dispersion relations of magnetic excitations have been 
determined by inelastic neutron scattering experiments. ^'^'^ 
The excitation spectrum consists of two distinct contributions: 
the dispersive modes at low energies below about 10 meV and 
the high energy mode around 24 meV with only weak dis- 
persion. These modes are assigned respectively to spin wave 
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Fig. 1. (Color online) The crystal structure of Cu2Fe2Ge40i3 projected 
onto the ab plane. The unit cell is shown by a solid line. 
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excitations of the Fe chains and singlet-triplet excitations of 
the Cu dimers. This assignment has been supported by a re- 
cent theoretical calculation of the magnetic excitation spec- 
trum.^- A set of exchange parameters has also been evalu- 
ated from analyses of the dispersion relations. The dominant 
interactions are all antiferromagnetic consistent with the mag- 
netic structure and are obtained as follows: the intradimer Cu- 
Cu interaction 7cu = 24 meV, the intrachain Fe-Fe interaction 
/pe = 1-60 meV, the interchain Fe-Fe interaction J^^ = 0.12 
meV in the c direction, and the Cu-Fe interaction 7cu-Fe = 
2.54 meV. 

One of the peculiarities of the magnetic excitations in 

Cu2Fe2Ge40i3 is that the Cu dimers and the Fe chains con- 
tribute almost independently to the excitation spectrum at dif- 
ferent energies, despite the fact that the coupling Jcn-ve be- 
tween the Cu and Fe subsystems is comparable to the intra- 
chain coupling /pe- Eventually, the low energy modes may 
consistently be explained by a spin wave theory for the Fe 
chains with the effective Fe-Fe coupling 7eff = 0.09 meV via 
the Cu dimer as well as the interchain coupling J'^^ in the c 
direction..^ Masuda et al. have pointed out that such a separa- 
tion of energy scales of the magnetic excitations between the 
two subsystems arises from a hierarchy of the exchange inter- 
actions 7cu » -^Fe. -^cu-Fe-^'^ The Suppression of the effective 
Fe-Fe interaction 7eff is interpreted as resulting from a non- 
magnetic nature of the Cu dimers which have an excitation 
energy much larger than the relevant energy scales of the Fe 
chains. 

Although the magnetic excitation spectrum of 
Cu2Fe2Ge40i3 may be understood as being contributed 
by effectively-decoupled Cu dimers and Fe chains, it is 
highly desirable to identify the modes more definitely using a 
microscopic probe complementary to the neutron scattering. 
Nuclear magnetic resonance (NMR) and nuclear quadrupole 
resonance (NQR) are best suited for such a purpose, be- 
cause they can probe local spin state and excitations quite 
sensitively. In this paper, we report on the results of NMR 
and NQR measurements of ^^Cu and ^^Cu nuclei (both 
having nuclear spin / = 3/2) in Cu2Fe2Ge40i3 and in an 
isostructural compound Cu2Sc2Ge40i3 in which Fe chains 
are all replaced by nonmagnetic Sc chains.'' Cu2Sc2Ge40i3 
contains Cu dimers which are magnetically well isolated 
from each other and have a spin gap of about 24 meV.'^ '"^ 
The room- temperature lattice parameters are a = 12.336 A, 
= 8.703 A, c = 4.888 A, and p = 95.74°." Our measurements 
of the ^^Cu nuclear spin-lattice relaxation rate in the two 
compounds reveal that Cu^^"^ and Fe^"^ spin fluctuations in 
Cu2Fe2Ge40i3 are decoupled to contribute additionally to 
the relaxation rate at the Cu site. The dynamics of Cu^^"^ spins 
in Cu2Fe2Ge40i3 is characterized by a spin gap identical to 
that in Cu2Sc2Ge40i3, which allows us to unambiguously 
identify the high energy mode as magnetic excitations of the 
Cu dimers. We also find from the relation between directions 
of the ordered moment of Cu and the intemal field at the 
Cu nuclear site that the electronic orbital of a Cu^^"^ ion in 
Cu2Fe2Ge40i3 has predominantly a d(x^ - y^) character. 

In §2 we describe the experimental procedures. The results 
are presented and analyzed in §3, where the characteristics of 
Cu^"^ and Fe^"^ spin fluctuations in Cu2Fe2Ge40i3 are deduced 
from the nuclear spin-lattice relaxation rate of Cu. The inter- 
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Fig. 2. Cu NQR spectra in Cu2Fe2Ge40n (solid line) and Cu2Sc2Ge40i3 
(dasiied line) at 50 K. Intensities are normalized at the pealc values. 



site are also determined in §3 through an analysis of the Cu 
NMR spectrum in the antiferromagnetic state. In §4 we dis- 
cuss an electronic state of the Cu^^ ion in the Fe compound 
and the nuclear spin-lattice relaxation rate at high tempera- 
tures in both the Fe and Sc compounds. Summary of the paper 
will be given in the final section 5. 

2. Experiments 

Single crystals of Cu2Fe2Ge40i3 were grown by the float- 
ing zone method.^ The dimension of the crystal used in the 
experiments was 5 x 3.5 x 3.2 mm-'. A polycrystalline sam- 
ple prepared by the solid-state reaction technique" was used 
for measurements in Cu2Sc2Ge40i3. NQR and NMR exper- 
iments were performed with a phase-coherent type pulsed 
spectrometer at zero external magnetic field. Standard n/l-n 
two-pulse sequence was used to excite spin-echo signals. The 
Cu NQR spectra were taken by Fourier-transforming the spin- 
echo signal with the step-sum technique.'^ The nuclear spin- 
lattice relaxation rate was measured by the inversion recovery 
method. The Cu NMR spectrum in the antiferromagnetic state 
of Cu2Fe2Ge40i3 was taken by recording an integrated inten- 
sity of the spin-echo signal point by point in the frequency 
range from 3.6 to 100 MHz. 

3. Results and Analyses 

3.1 Cu NQR spectrum 

A pair of sharp resonance lines was observed at zero field in 
the paramagnetic state of Cu2Fe2Ge40i3 (Fig. 2). We identi- 
fied the signal as Cu NQR from a unique crystallographic site 
because a ratio between the peak frequencies agrees with that 
between the electric quadrupole moments of *^Cu and ^^Cu, 
63g/65g ^ 2 081. This is consistent with flie crystal struc- 
ture in which aU the Cu atoms occupy equivalent positions.^ 
The ^^Cu NQR frequency determined from the peak position 
of the resonance line is 34.90 MHz at 50 K. The Unes are 
broadened by inhomogeneous distribution of the electric field 
gradient (EFG) at the Cu site. The full width at half maximum 
(FWHM) is 149 kHz at 50 K for ^^Cu. 

In Fig. 2 the Cu NQR spectrum in Cu2Sc2Ge40i3 is also 
shown for comparison. The ^^Cu NQR frequency is 34.91 
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Fig. 3. Temperature dependence of the *^Cu NQR frequency in 
Cu2Fe2Ge40i3 (solid circles) and Cu2Sc2Ge40i3 (open squares). 



same temperature within experimental accuracies. On the 
other hand, the hne is broader in Cu2Sc2Ge40i3 with the 
FWHM of 206 kHz at 50 K for ^Scu. This indicates that 
the local structure around Cu is more homogeneous in the 
Cu2Fe2Ge40i3 crystal specimen. 

The temperature dependence of the ^^Cu NQR frequency 
*^vnqr in the two compounds is almost identical as shown in 
Fig. 3, although the data for Cu2Sc2Ge40i3 are limited up to 
90 K due to fast nuclear spin relaxation processes, ^^vnqr in 
Cu2Sc2Ge40i3 is independent of temperature below about 50 
K. As the temperature is raised, ^'^vnqr continues to increase 
up to 260 K in Cu2Fe2Ge40i3. This is opposite to a tendency 
of the EFG to decrease with increasing temperature due to 
lattice thermal expansion. Such a behavior of the Cu NQR 
frequencies has been reported in some Cu oxides including 
high-Tc superconductors, and is attributed to a dominant con- 
tribution of on-site electronic orbitals to the EFG which is 
different in sign from the contribution of surrounding ions. 

3.2 Nuclear spin-lattice relaxation 

The nuclear spin-lattice relaxation rate 1/7'i,nqr at the 
Cu site in Cu2Fe2Ge40i3 was determined as a time con- 
stant of exponential recovery of the nuclear magnetization. 
Non-exponential recovery was observed in Cu2Sc2Ge40i3 
below 45 K due to contributions of paramagnetic impurities. 
1/^1, NQR was then determined by fitting the nuclear magneti- 
zation M{t) to the stretched exponential form which incorpo- 
rates both the impurity and intrinsic relaxations, 



M(f) = M^\\- po exp (- ^|t/Tc - ^/Tlnqr)], (1) 

where Moo is the magnetization in thermal equilibrium, po 
is a parameter describing a degree of inversion, and 1 /tc is 
the relaxation rate induced by paramagnetic impurities. The 
impurity relaxation rate I/t^ is almost temperature indepen- 
dent and takes a value of about 1 s ' for ^^Cu. In both com- 
pounds the isotopic ratio of I/T^i^nqr between *^Cu and ^^Cu 
agrees with the ratio of the square of the gyromagnetic ra- 
tios (^3y/65^-)2 Yyifjjjjj experimental accuracies. The relaxation 
process is thus magnetic in origin. 

Fi gure 4 shows the temperature dependence of 1/7^i,nqr 
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Fig. 4. Temperature dependences of the nuclear spin-lattice relaxation 
rate l/ri_NQR at the *^Cu site in Cu2Fe2Ge40i3 (soUd circles) and 
Cu2Sc2Ge40i3 (open squares). Dashed line stands for Tn of Cu2Fe2Ge40i3. 
The inset is a plot of l/ri_NQR(Fe) versus 1 /7'i_nqr(Sc) with temperature the 
implicit parameter. SoUd Une is a fit of l/7'ij,jQR(Fe) above 50 K. 



1 /T'i.nqrCFc) in Cu2Fe2Ge40i3 increases gradually with in- 
creasing temperature well above Tn, exhibiting a tendency to 
saturate at high temperatures. On the other hand, it increases 
divergently on approaching T^. This signals critical slowing 
down of electronic spin fluctuations toward long-range mag- 
netic order. 

As a clear contrast, a monotonous and rapid decrease of the 
relaxation rate I/Ti nqrCSc) in Cu2Sc2Ge40i3 was observed 
at low temperatures. This reflects a nonmagnetic ground state 
with a finite spin gap as expected for spin dimers. Fig- 
ure 5 shows l/ri,NQR(Sc) as a function of inverse tempera- 
ture 1/r which evidences a thermally-activated behavior of 
l/ri,NQR(Sc). By fitting the data below 50 K to the formula 
I/^i.nqrCSc) = cexp(-As/r), we determined the constant 
c and the spin gap As to be (2.3 ± 0.1) x 10^ s"' and As = 
269 ± 10 K, respectively. The value of As agrees well with 
those estimated from susceptibiUty, '^^Sc-NMR, and neutron 
scattering measurements.^^' 

It is instructive to compare the temperature dependences 
of l/ri,NQR(Fe) and l/ri,NQR(Sc) to get insight into how the 
low-energy magnetic excitations of Cu dimers are affected 
by Fe chains in Cu2Fe2Ge40i3. The inset of Fig. 4 is a plot 
of l/ri,NQR(Fe) versus l/ri_NQR(Sc) with temperature the 
imphcit parameter. It is clear that l/ri,NQR(Fe) scales with 
1/7^i,nqr(Sc) except at temperatures very close to Tn. Above 
50 K, l/ri,NQR(Fe) is expressed as 

^- = fl H- b ^- , (2) 

ri,NQR(Fe) ri,NQR(Sc) 

where a and b are temperature-independent constants. The 
values of a and b were determined by the least-square fitting 
to be (2.67 + 0.05) x 10^ s'^ and 0.213 ± 0.004, respectively. 
Deviation from the scaling near is due to critical slowing 
down. 

Equation (2) demonstrates two important characteristics of 
l/7'i,NQR(Fe): first, it has a temperature-independent contri- 
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Fig. 5. Nuclear spin-lattice relaxation rate l/Tj nqr(M) at the *'Cu site in 
Cu2Sc2Ge40i3 (M = Sc; solid squares) and Cu2Fe2Ge40i3 (M = Fe; dotted 
line) plotted against l/T. Solid line is a fit of l/Ti nqr(Sc) to the activation 
law with a spin gap As = 269 K. Solid triangles are the contribution of Cu^+ 
spins (l/Ti NQR(Fe))ci, to l/Ti NQR(Fe). Dashed line is drawn as a guide to 
the eyes to demonstrate an activated behavior of (l/Ti NQR(Fe))cu with the 
same spin gap as that of l/Ti nqr(Sc). 



bution which is absent in Cu2Sc2Ge40i3, and second, the 
temperature-dependent part is qualitatively the same as the 
temperature dependence of 1 /Ti nqr(Sc) except near Tn. The 
temperature-independent term a should be interpreted as a 
contribution of Fe"'^ spin fluctuations to l/ri,NQR(Fe), be- 
cause it is the term that remains as 1 /Ti nqrCSc) 0. The 
second term of eq. (2) is regarded as the relaxation rate result- 
ing from Cu^^ spin fluctuations under the influence of the Fe 
chains. 

It is worth noting that the contribution of Fe^^ spins is 
additional to that of Cu^^ spins. This means that the Fe^^ 
and Cu^^ spins contribute independently to l/ri_NQR(Fe). The 
Fe^^ spins affect the contribution of Cu^+ spins only through 
a numerical factor b in eq. (2), without changing a thermally- 
activated character of Cu^^ spin fluctuations. In Fig. 5 we 
plotted the contribution of Cu^+ spins to l/Ti NQR(Fe) eval- 
uated as 



1 



1 



- a (3) 

v7'i,NQR(Fe);py ri,NQR(Fe) 

against l/T. Raw values of 1 /ri,NQR(Fe) were also plotted for 
comparison. Clearly, (1 /Ti NQR(Fe))cu parallels 1 /Ti nqr(Sc) 
above about 50 K. This indicates that the temperature depen- 
dence of (1 /ri_NQR(Fe))cu is characterized by a spin gap iden- 
tical to that of 1 /ri,NQR(Sc). 

All these observations are consistent with the results of 
the neutron scattering experiments.^- ^-^ In Cu2Sc2Ge40i3, 
a singlet-triplet excitation mode of Cu dimers is observed 
at high energies around 24 meV corresponding to the in- 
tradimer exchange interaction. This mode is persistent in 
Cu2Fe2Ge40i3 at the same energy and is associated with Cu 
dimers. Our result of (l/T) NOR(Fe))ru confirms this assign- 
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Fig. 6. (Color onHne) Cu NMR spectmm in the antiferromagnetic state of 
Cu2Fe2Ge40i3 at 4.2 K and T. Solid line is a simulation of the spectrum 
based on the parameters given in the text. Dashed and dotted lines are the 
contiibutions of *'Cu and *^Cu to the spectrum, respectively. Gaussian distri- 
bution of the internal field with the FWHM of 3. 1 MHz for '^Cu was assumed 
to simulate the spectrum. 



ment. The gapped excitations of Cu dimers is thus left un- 
changed by Fe chains. On the other hand, excitation modes at 
low energies below about 10 meV are viewed as spin waves 
propagating only within the Fe subsystem. The low-energy 
modes contribute to the excitation spectrum independently of 
the dimer mode as if the Cu dimers have nothing to do with the 
dynamics of the Fe chains. The fact that the Fe^*^ spins con- 
tribute independently of the Cu^+ spins to 1 /ri,NQR(Fe) must 
be related to this feature of the excitation spectrum which 
would result from a large difference of relevant energy scales 
between the two subsystems. 

While the temperature dependence is governed by a spin 
gap which is common to both compounds, the absolute val- 
ues of l/ri,NQR(Sc) and (l/ri,NQR(Fe))cu are scaled by the 
numerical factor b. This factor is determined primarily by a 
ratio of the hyperfine coupling constant at the Cu site between 
the two compounds. Likewise, the value of a contains various 
information such as the dynamics of Fe^^ spins and the hyper- 
fine interaction between the Cu and Fe sites. We will return to 
these subjects in the next section. 

3.3 Cu NMR spectrum in the antiferromagnetic state 

We observed zero-field resonance of Cu nuclei in the an- 
tiferromagnetic state of Cu2Fe2Ge40i3. The lines are largely 
shifted from the NQR frequencies in the paramagnetic state, 
indicating appearance of internal magnetic field at the Cu site 
due to ordering of Cu and Fe magnetic moments. 

The spectrum is shown in Fig. 6. At 4.2 K and T, three 
groups of resonance lines, cc' and hh', are observed in the 
frequency range from 3.6 to 80 MHz. Such a spectrum may 
be interpreted as a superposition of the NMR lines of the two 
isotopes ^^Cu and ^^Cu split by the electric quadrupole inter- 
action. Indeed, it is confirmed from the following analysis that 
the lines c and c' correspond respectively to the central tran- 
sitions of ^^Cu and ^^Cu. whereas the lines / and h (/' and h') 
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Table L Observed peak frequencies and assignment of tlie resonance lines 
of the Cu NMR spectmm. Frequencies calculated using a set of parameters 
given in the text are also shown for comparison. Calculated intensities are 
normalized to the ^'Cu central line c. 



line 


Vobs (MHz) 


Veal (MHz) 


/cal 


/ 


5.8 ±0.2 


5.67 


0.768 


*3Cu c 


37.8 ± 0.2 


37.68 


1 


d 




43.35 


0.162 


h 


72.5 ± 0.3 


72.40 


0.866 


r 


9.8 ±0.3 


10.05 


0.430 


•^^Cu c' 


40.8 ± 0.2 


40.78 


0.571 


d' 




50.82 


0.020 


h' 


72.5 ± 0.3 


72.59 


0.443 



are the quadrupole satellites of ''^Cu (*^Cu). The overlap of 
the high-frequency satellites h and h' is the reason for the re- 
duced number (five rather than six) of the resolved peaks. The 
peak frequencies and the assignment of the lines are summa- 
rized in Table I. 

NMR frequencies in the antiferromagnetic state under zero 
external field are generally determined by the Hamiltonian 
consisting of both the Zeeman and quadrupole interaction 
terms "TYz and "TYq: 



'Hz = -jni ■ Bint 



h COS 6* -H - sin e [l+e-'^ + /_e'<*) 



(5) 



hVr 



Vq = 



6 

"ieQVzz 



I2 



with 



2/(2/- 1)/! 

Here we refer to the principal frame of the EFG tensor y is 
the gyromagnetic ratio of a nucleus, Bint is the internal field 
at the nuclear site, Q is the nuclear quadrupole moment, Vzz 
is the Z component of the EFG tensor at the nucleus, and rf is 
the asymmetry parameter of the EFG defined as t; = iVxx - 
Vyy)IVzz with V^^ = ^'^V|^^? in = X,Y,Z). We take \Vxx\ 
< \Vyy\ ^ I Vzz I by convention so that Q < j] < 1. Polar and 
azimuth angles 9 and (p specify the orientation of Bint with 
respect to the principal axes of the EFG tensor 

In order to determine the five parameters Bint, S, (p, vq 
and which best reproduce the observed resonance frequen- 
cies, we numerically diagonalized the Hamiltonian eq. (4) 
and calculated the resonance frequencies by adjusting param- 
eter values. This is because in the present case the Zeeman 
and quadrupole interactions are comparable in magnitude, so 
that the perturbation theory cannot give resonance frequencies 
with desired accuracy."* 

In searching for the parameter values, we used the value 
s^vnqr - 34.90 MHz at 50 K as a value at 4.2 K, because the 
NQR frequency vnqr in Cu2Fe2Ge40i3 is considered to be 
temperature independent below 50 K judging from the tem- 
perature dependence of ^^vnqr in Cu2Sc2Ge40i3 shown in 
Fig. 3. Since vnqr for nuclei with / - 1^12 is given as 



Vnqr = V ^ "3"' 



(7) 
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(4) Fig. 7. (Color online) (a) Comparison of the observed quadrupole split- 
tings Aotis (dotted lines) with vq calculated as a function of from ^^vnqr 
= 34.90 MHz. Solid and dashed lines represent ^^vq and ^''yq, respectively. 
The condition A^tis < yq is required for both isotopes, (b) Comparison of the 
observed quadrupole splitting Aob, (dotted lines) with the splitting A calcu- 
lated as a function of d for = (solid lines), 0.30 (dashed lines) and 0.46 
(dotted dashed lines) using eqs. (7) and (8) with '^vnqr = 34.90 MHz and 
= 0°. For a given value of the solution of A(6) = A^bs yields a polar 
angle 6 of the internal field. 

(6) 



63, 



Vq can be determined from vnqr by setting 77. vg is then 



set by Vq using the isotopic ratio of the quadrupole moments 
as '''^vq - ^^vq{^^QI''^Q). Values to be determined then re- 
duce to four; Bint, S, and 77 which are common to both iso- 
topes. 

Although the perturbation theory cannot be adopted to ob- 
tain accurate values of Bint, <P and t], the results might be 
used to find appropriate ranges of parameter values. Accord- 
ing to the second order perturbation theory, the quadrupole 
splitting A between the high- and low-frequency satelUtes of 
1-3/2 nuclei is given by'^ 



Vh - Vl 



— (3 cos^ ■ 
2 



1 +77 sin^ 61 cos 20), (8) 



where Vh and V] are resonance frequencies of the high- and 
low-frequency satellites, respectively. A satisfies a relation 

1+'? 



-Vq < A < Vq. 



(9) 



The upper and lower bounds of A correspond respectively 
to e = 0° and - (p - 90°. Therefore, |A| cannot exceed 
Vq = vnqr(1 + 77^/3)"'^^. Figure 7(a) compares the observed 
quadrupole spHttings "Aobs = 33.35 MHz and ''^Aobs =31.35 
MHz with " Vq and ^^vq calculated from " vnqr = 34.90 MHz 
as a function of 77 using eq. (7). The conditions ^^Aobs < ^^vq 
and *^Aobs < ^^vq are satisfied simultaneously when 77 < 0.46 



which we resard as an aDnroDriate ranee of 77. 
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Next we consider a possible range of 0. Figures 7(b) shows 
plots of A versus 6 for several values of rj calculated using 
eqs. (7) and (8) with "vnqr = 34.90 MHz and (P = 0°. In 
this plot, the solution of A(ff) - Aobs yields a polar angle 
of the internal field for a given value of t]. For example, the 
calculated splitting A for ''^Cu crosses with *^Aobs at - 9° 
when 77 = 0.30, yielding a polar angle of 9° for this value of 
77. The observed splitting ^^Aobs (''^Aobs) hence corresponds to 
a value of 61 in a range 6° < < 1 1° (0° < 6* < 9°) depending 
on T]}'^ The case of 0° (not shown) gives qualitatively the 
same results. The range 0° < < 11° would therefore be suf- 
ficient for the parameter search. Note that this range of has 
an important implication on the direction of Bjnt: the internal 
field is nearly parallel to the Z principal axis of the EFG ten- 
sor, being tilted by at most ~10°. A possibility of ~ 90° can 
safely be excluded because for rj < 0.46 the splitting is at most 
(1 + t])vq/2 = 0.73vq = 24.62 MHz for ''^Cu which is much 
smaller than ^"'Aobs- 

The magnitude Z?int of the internal field may roughly be 
evaluated from the frequency of the central transition, neglect- 
ing the second order effect of the quadrupole interaction. The 
resonance frequency Vc is then given as Vc = yBint- Using the 
values Vc = 37.8 MHz (40.8 MHz) and y = 11.285 MHz/T 
(12.089 MHz/T) for ^^Cu (''^Cu), we obtain Bim = 3.35 T 
(3.37 T) as an approximate value. 

Having had appropriate ranges of the four parameters, we 
searched the values of Bint, &, <P which minimize a quan- 
tity = ^/[(nobs - v/,cai)^/o-^] by solving numerically the 
Hamiltonian eq. (4). Here v;,obs (>'/,cai) is the observed (calcu- 
lated) resonance frequency of the i-th line, cr, is the experi- 
mental uncertainty of v,;obs, and the sum runs over all the ob- 
served lines /, r, c, c', h and h' listed in Table I. Parameter 
ranges investigated were 3.1 < Bim < 3.6 T, 0° < 6* < 11°, 0° 
< < 90° and < /; < 0.46.^' We found that the resonance 
frequencies are insensitive to for such a small value of 0. We 
therefore optimized Bint, and rf by fixing (p to some values. 

We successfully reproduced the observed resonance fre- 
quencies by properly setting Bint, 0, (p and 77 within the above 
ranges of parameters. The optimized values of and rj depend 
slightly on 0: they change from 4.1° and 0.283 for <p = 0°, to 
3.4° and 0.296 for 4>^90°. The optimized value of Bint, on the 
other hand, stays constant to take a value of 3.39 T when is 
varied from 0° to 90°. These sets of parameters give almost 
the same resonance frequencies and the minimum value - 
1.7. The resulting v,,cai's are listed in Table I. 

We now wish to know what value of </> best describes 
the observed NMR spectrum. For that purpose we calcu- 
lated the intensity of each resonance line using the eigen- 
states of the Hamiltonian eq. (4) for the optimized sets of 
parameters. Intensity 1^^ of the resonance line arising from 
transition between the states \a} and is given as lajs 
Ny^fi^\{a\I ■ Bi|y6)p, where is the abundance of a nuclear 
species and Bi is an rf exciting field. We take a powder aver- 
age of I„js because the direction of Bi with respect to the EFG 
principal axes is unknown. 

Figure 8 shows NMR spectra calculated using the opti- 
mized values of Bint, & and rj for (p - 0° and those for <p = 90° 
presented above. The most striking feature of the calculated 
spectra is existence of weak resonance lines d and d' not ob- 
served experimentally. These lines correspond respectively to 
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Fig. 8. (Color online) Cu NMR spectra calculated using the sets of param- 
eters (a) B,n, = 3.39 T,7j = 0.283, 9 = 4.1° and (p = 0°, and (b) Bi,,, = 3.39 T, r/ 
= 0.296, e = 3.4° and = 90°. Solid (dashed) lines labeled by letters /, c and 
h {l',c' and h') represent low-frequency satellite, central, and high-frequency 
satellite lines of *'Cu (*^Cu), respectively. Lines d and d' correspond respec- 
tively to transition of '^Cu and *'Cu which is forbidden in the limit of strong 
Zeeman interaction. Intensities are normalized to the central line c. 



in the limit ITYzl ^ I'KqI and become more intense as <p varies 
from 0° to 90°.^^ Of the two, the line d seems strong enough 
to be detected by experiments. The most probable explanation 
for the lack of observation is that the line d is masked by more 
intense lines such as c and c'. 

The line d would then appear, if observed, as an asymmet- 
rical tail in the high frequency side of the central lines c and 
c' which becomes more pronounced as <p increases. Since we 
have no definite sign of such a structure in the observed spec- 
trum, a set of parameters for = 0° which yields the minimum 
intensity of the lines d and d' seems best reproduce the exper- 
imental result. The fact that the line I is more intense than the 
line c when <p - 90° favors the (p - 0° parameter set as well, 
because this is opposite to the observation. We finally get the 
following set of parameters: 

Bint = 3.39 + 0.01 T, 

61 = 4.1° +0.1°, 



63 



■■ 0°, (10) 

vg = 34.44 + 0.04 MHz, 

77 = 0.283 + 0.005. 

Errors in Bint, O and 77 were evaluated as a range of parameters 
giving < 1.8. The obtained value of Bint is modestly large 
for 3d transition metal ions and may be attributed to reduction 
of the ordered moment of Cu due to quantum fluctuations. 

The intensities of the resonance lines calculated using the 
above parameters are given in the last column of Table I. We 
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Gauss functions based on the resonance frequencies and in- 
tensities listed in Table I. The result is shown in Fig. 6 for 
comparison with the experiments. 

4. Discussion 

4.1 Internal field and 3d electronic orbital 

In the previous section, we determined the direction as 
well as the magnitude of the internal field at the Cu site in 
Cu2Fe2Ge40i3. The internal field at the nucleus belong- 
ing to a magnetic ion comes mostly from the ordered moment 
m of that ion via the (on-site) hyperfine coupling as 

Am. 



B„ 



(11) 



A is the hyperfine tensor which is determined by an electronic 
state of the magnetic ion. Therefore, the internal field carries 
information not only on the ordered moment but also on the 
electronic orbital. In this subsection, we will discuss the in- 
ternal field and the 3t/ electronic state at the Cu site in detail, 
focusing on the relative orientation of the internal field and 
the ordered moment. 

To proceed further, we must specify the orientation of the Z 
principal axis of the EFG tensor with respect to the crystalline 
axes which is as yet undetermined. We assume here that the Z 
axis is normal to the CU2O2 dimer plane based on the fourfold 
planar coordination of Cu, although there is some distortion 
from the ideal square planar geometry. This axis is also taken 
as the Z principal axis of the hyperfine tensor. The angle be- 
tween the Z axis and the ordered moment nicu of the Cu atom 
is then calculated to be 60° (or 120° for another sublattice) 
using the result of the neutron diffraction experiments.^ This 
is a clear contrast to our observation that the internal field Bint 
at the Cu site is canted scarcely from the Z axis with = 4.1°. 
Bint and iricu are therefore noncollinear, making a large angle 
of about 60° or more between them. 

It is clear from eq. ( 11 ) that Bint and m become noncollinear 
when A is anisotropic and m is parallel to none of the principal 
axes of A. In the present case, the hyperfine interaction of the 
Cu^^ ion is highly anisotropic because a 3d hole occupies a 
spatially anisotropic orbital such as dix^ - y^) and £/(3z^ - r^). 
Regarding the ordered moment iricu, it is not parallel to any 
of the principal axes. The reason why mcu and Bint at the Cu 
site are noncollinear is understood in this way. 

We now ask whether the anisotropic hyperfine interaction 
characteristic of the Cu^^ ion is able to account for the small 
canting angle of Bint. The most stringent test would be to cal- 
culate Bint from the experimentally-determined A and iricu- 
Unfortunately, A at the Cu site in Cu2Fe2Ge40i3 is undeter- 
mined because of difficulty in detecting NMR signals in the 
paramagnetic state. Here we present a model calculation of 
the dependence of Bint on the direction of iricu using A typical 
for the Cu"^ ion. 

In calculating Bint we assume for simplicity that the Cu site 
has tetragonal symmetry about the Z axis. The ground state 
orbital is then either d(x^ - y^) or d{3z^ - r^). In reality, the 
point symmetry at the Cu site is lower than tetragonal, which 
will cause mixing of orbitals with various symmetry and re- 
sultant deviation of the principal axes from their ideal orien- 
tations. Nevertheless, our simplified model well describes the 
observed features of Bint if we take d(x^ -y^) as a ground state 
orbital. 
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Fig. 9. (Color online) Relation between the angles 9b and 0,,, for d(x^ -y^) 
(solid line) and d(3z^ - r^) (dashed line) orbitals. The inset shows Bi„t as a 
function of (9,^ for the two orbitals with mcu = 0.38 /jg. Dotted lines in the 
main panel and the inset correspond to 6^ of 60° and 120° . 



sufficient to work in the plane containing mcu and the Z axis. 
We take the X' axis in this plane perpendicular to the Z axis. 
Defining 0^ (6b) as an angle between iricu (Bint) and the Z 
axis, we have 



mcu = (mcu sin 0m, mcu cos 0,n), 
Bint = (Bint sin 6Ib, Bint cos %), 



(12) 
(13) 



where mcu - Inicul- Utilizing eq. (11), we obtain the following 
relations: 

Bint = mcJAl cos2 0,, + Al sin^ 0^)"\ (14) 

tan^B = — tanffm- (15) 
^11 

Here Ay - Azz and Aj_ = Ax'x' are the principal values of 
A. They are estimated as Ay - -24.1 T/^ub and Aj_ - 2.71 
T/;Ub for d(x^ - A|| = 10.2 T/^ub and A^ = -12.2 T/fi^ for 
d(3z^ - r^).P Using these values, we calculated Bint and ^b 
as a function of 9^ for d{x^ - y^) and d{3z?' - r^). 

Figure 9 shows the angle 0^ plotted against 6^. Obviously, 
we always have a negative slope of 9^ with respect to 9^, i.e., 
d9B/d9^ < for both the orbitals. In other words. Bint rotates 
counterclockwise when iricu rotates clockwise as illustrated in 
Fig. 10. This comes from the fact that Ay and A± are different 
in sign. The dependence of 9^ on 9^ is, however, qualitatively 
very different between the two orbitals. 

We examine the result for d(x^ - y^) first. When iricu is di- 
rected along the Z axis, i.e., 9„^ - 0°, Bint also lies along the 
Z direction, being antiparallel to iricu with 9b - 180° (Fig. 
10(a)). This reflects the fact that A|| < for d(x^ - /). The 
angle 9^ decreases gradually with increasing 9^, yet the resul- 
tant canting of Bint from the Z direction remains small, being 
only within 11° (9b > 169°) for 9^ up to 60°. Significant 
change of 9b occurs in a range 80° < 0m ^ 100° where more 
than 110° change of 9b takes place due to only 20° change 
of 9in- nicu and Bint become collinear again along the X' axis 
with ^m ^9b^ 90° (Fig. 10(c)). 
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(b) 9^ = 60° 
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Fig. 10. (Color online) Change of the direction of the internal field at the 
Cu site accompanied by a rotation of the Cu magnetic moment. Solid arrow 
represents the Cu magnetic moment. Dashed and dotted arrows are the inter- 
nal fields for d{x^ - y^) and d(3z~ - r^) orbitals, respectively. The angle 
between the magnetic moment and the Z axis is (a) 9°, (b) 69°, (c) 99°, (d) 
129°, and(e) 189°. 



with 0^. Since Ay > 0, iricu and Bjnt are parallel to each other 
when din = 0°. Ob varies almost linearly with increasing 6^, 
satisfying the approximate relation ~ -dm- Along the X' 
axis. Bint is oriented antiparallel to iricu- 

When 6m - 60° or 120° which is the case in Cu2Fe2Ge40i3, 
we have 0b - 169° or 11° for d{x^ - It should be noted 
here that there is no distinction between the two sets of the 
angles (0m, 0b) = (60°, 169°) and (120°, 11°): they are equiva- 
lent because of twofold symmetry about the X' axis. See Figs. 
10(b) and 10(d). In either case the deviation of Bjnt from the 
Z direction is only 11° despite the large canting angle (60°) of 
mcu- On the other hand. Bint for d{3z^ - r^) is directed away 
from the Z direction to have a large canting angle of 64° (0b 
= -64° or -116°). 

The dependence of Bint on 0m is also very contrasting be- 
tween the two orbitals as shown in the inset of Fig. 9. Bint for 
d{x^ -y^) depends strongly on 0m- for mcu - 0.38 jiB, it varies 
from 9 to 1 T as 0n, varies from 0° to 90°. In contrast, Z?int 
for d(3z^ - r^) changes only slightly with 0m, ranging from 
3.9 to 4.6 T. When 0m = 60° and 120°, two orbitals give acci- 
dentally the same value Bin, = 4.6 T. This is not very different 
from the observed value of 3.39 T. Notice also that the large 
canting angle of mcu as well as the reduced value of otcu is 
responsible for the modestly large value of Bint for d{x^ -y^)- 

It is apparent from the above discussion that anisotropy of 
the hyperfine interaction of the d{x^ - y^) orbital accounts 
qualitatively for the small canting angle of Bint which arises 
from mcu being tilted away from the symmetry (Z) axis. This 
suggests a dominant d{x^ - y^) character of the Cu^+ ions in 
Cu2Fe2Ge40i3. The d{x^ -y^) character would be understood 
by remembering the fourfold planar coordination of Cu. The 
fourfold, square planar coordination of Cu by O is derived 
from a CuOg octahedron by removing the apical oxygens. The 
point symmetry at the Cu site is tetragonal but the d{x^ - y^) 
orbital is energetically favored because there is no apical oxy- 
gen. Note also that this orbital extends to the directions of 
Cu-O bonding, allowing the Cu^+ ion to exchange interact via 



common edge of the CUO4 squares. The d{x^ - y^) nature 
of the ground state orbital seems thus compatible with the 
strong intradimer interactions observed in Cu2Fe2Ge40i3 and 

CU2SC2Ge40i3. 

There are some quantitative discrepancies between the ex- 
perimental results and the model calculations, especially in 
the magnitude Bint of the internal field. This may be attributed 
to mixing of the d(3z^ - r^) orbital to the ground state allowed 
by distortion from the ideal square planar geometry which is 
manifested by a moderate asymmetry 77 = 0.283 of the EFG 
tensor A preliminary model calculation shows that 10 to 20 
% mixing of the c/(3z^ - r^) orbital improves dramatically 
the quantitative discrepancies, yielding a canting angle of Bint 
within 6° from the Z axis and Bint of 3.3 to 3.9 T. 

4.2 Nuclear spin-lattice relaxation at high temperatures 

Nuclear spin-lattice relaxation rate is one of the most im- 
portant quantities measured by NMR and NQR, because it re- 
flects low-energy magnetic excitations quite sensitively. Here 
we discuss the nuclear spin-lattice relaxation at the Cu site at 
high temperatures in both Cu2Sc2Ge40i3 and Cu2Fe2Ge40i3. 

First we discuss the relaxation rate in Cu2Sc2Ge40i3. In the 
system with localized magnetic moments, the nuclear spin- 
lattice relaxation rate of / = 3/2 nuclei measured by NQR is 

24 25 

given as • 

1 y^SiS + 1) 



rioo,NQR V 2 3(3 + ?7-)Wex 

X [(3 + iifAl^ + (3 - iTfA\y + Ait'aI^] (16) 

in the high-temperature limit. Here S is the electronic spin, 
A^^'s (p - X, Y,Z) are the principal values of the hyperfine 
tensor Wex is the exchange frequency which is given for an 
isolated spin dimer as 



(17) 



where Jo is the intradimer exchange interaction. If we assume 
axial symmetry of the hyperfine tensor Axx - Ayr = A^ and 
Azz - A|| for simplicity, eq. (16) reduces to 

[1^2y^AlS(S + 1) 
2' 



1 



loo, NQR 



/('?,'-) = 3 + 



2T]\r^ - 1) 



(18) 



liifli din 



3 +T]^ 

Here r = Ay/Aj^, and the function /(tj, r) describes effects of 
asymmetry 77 of the EFG and the anisotropy r of the hyperfine 
coupling. Putting S = 1/2, - 24 meV into eq. (17), we ob- 
tain Wex = 2.6 X 10'^ s"'. If we take the d(x^ -y^) orbital state 
for a Cu^^ ion, we have r = -8.89. We then obtain I/Tioo.nqr 
= 2.5 X lO'* S-' for ^^Cu from eq. (18) using 77 = 0.283. Note 
that for T] = 0.283 and r = -8.89, f(i],r) = 7.06 which is 
much larger than the value of the function for cases with 77 
= and/or r - 1. This value of l/rioo,NQR seems to be in the 
right orders of magnitude because 1 /ri_NQR(Sc) is expected to 
reach around a value of 10^ s"' in the high-temperature limit 
judging from its temperature dependence shown in Fig. 5. 

In Cu2Fe2Ge40i3, we expect a contribution of the neigh- 
boring Fe^+ spin to 1 /Tiocnqr at the Cu site via the transferred 
hyperfine interaction in addition to the on-site one represented 
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Fig. 11. Local arrangement of magnetic atoms around Cu and definition of 
tiie exchange interactions Jo, Ji and J2. 



by eq. (16) or (18). I/Tioo.nqr at the Cu site in Cu2Fe2Ge40i3 
is hence expressed as 



1 



1 



1 



(19) 



7"loo,NQR \7"ioo,NQr/cu \7"loo,NQR/Fe 

where the first and second terms are the contributions of Cu^^ 
and Fe^"^ spins, respectively. Assuming the isotropic trans- 
ferred hyperfine interaction with a coupling constant B from 
the nearest neighbor Fe^^ ion, we have 

[l2y^AlScn(Scn + l)_ 



1 



^ri«,,NQR/cu 
/ 1 



3Wex,Cu 
„2d2( 



'-f{Tj,r), (20) 



2j'B'Sp,(Sp, + 1) 



3(0, 



(21) 

' loo,NQR/Fe 

Here 5cu = 1/2 (5 pe = 5/2) is the electronic spin of a Cu^"^ 
(Fe-'^) ion. Note that a factor 3 in the right hand side of eq. 
(21) comes from the fact that we are dealing with NQR. The 
exchange frequency o)ex,i of the electronic spin at the i-th site 
is generally given as 

\2 



= ^Z^/^.-i)(y)' (22) 



where S j is the electronic spin at the j-th site and Jjj is the ex- 
change interaction between the i-th and j-lh spins. Applying 
eq. (22) to the Cu^"^ and Fe^"^ ions in Cu2Fe2Ge40i3, we have 



"ex,Cu 



35 



(I) 



(23) 



(24) 



The exchange interactions Jq, Ji and J2 are defined as shown 
in Fig. 1 1 . We neglected the interchain Fe-Fe interaction (J^^) 
in the c direction because it is an order of magnitude smaller 
than Ji and J2. Using the values Jq - Jcu = 24 meV, Ji = 
Jpe = 1-60 meV, J2 = /cu-Fe = 2.54 meV determined by the 
inelastic neutron scattering experiments,^ we obtain a)ex,cu = 
2.7 X 10" S-' and Wex.Fe = 8.7 x lO'^ 

It is reasonable to consider the temperature-independent 
contribution of Fe^^ spins to l/ri NQR(Fe), i.e., the term a 
of eq. (2), as (l/rioo,NQR)Fe- By equating (l/rioo,NQR)Fe to a 
= 2.67 X 10^ s"\ the only unknown quantity in eq. (21), the 
absolute value of the transferred hyperfine coupHng B, is esti- 
mated as \B\ = 0.23 T/fiB. This is an order of magnitude larger 
than a classical dipole field from the surrounding Fe^"^ ions 
and is in a reasonable range for the transferred hyperfine in- 



The on-site contribution (1/7'ioo nqr)cu is estimated to be 
2.4 X 10'* s ' using Wex.Cu calculated above and assuming the 
dix^ -y^) orbital state. This is in good agreement with a value 
(l/7'ioo,NQR(Fe))cu ~ 2 X lO'* s"' evaluated by extrapolating 
the temperature dependence of (l/ri,NQR(Fe))cu to T — > 00 
utiUzing Fig. 5. 

Lastly we discuss the numerical factor b in eq. (2) relating 
l/ri,NQR(Sc) to (l/ri,NQR(Fe))cu- Since the temperature de- 
pendences of 1/7'i_nqr(Sc) and (l/7'i,NQR(Fe))cu are identical 
apart from the factor h, we may argue their values in the high- 
temperature limit using eqs. (18) and (20). The difference of 
the two quantities then comes primarily from the difference 
of the hyperfine coupling constant between the two com- 
pounds, because the exchange frequencies of the Cu^^ spin 
given by eqs. (17) and (23) are nearly the same. Utilizing 
the fact that 1 /Tioo^nqr A]_ and the obtained value of h, we 
evaluated a ratio of Aj_ between the Fe and Sc compounds as 
A_L(Fe)/Ax(Sc) ~ fci/^ ^ 0.46. 

It is known that a small and positive value of Aj_ for the 
dix^ - y^) orbital results from subtle balance of core polariza- 
tion, dipole and orbital fields which are different in sign and 
are an order of magnitude larger than the total A^}"^ On the 
one hand, we have in Cu2Fe2Ge40i3 the exchange interaction 
between Cu and Fe atoms which is absent in Cu2Sc2Ge40i3. 
It is thus probable that a microscopic process relevant to the 
Cu-Fe interaction modifies a wave function of the Cu atom in 
Cu2Fe2Ge40i3 from that in Cu2Sc2Ge40i3, altering the bal- 
ance of various contributions to A^. The factor of ~2 differ- 
ence of A J. between the two compounds might be attributed 
to such a change of the electronic state of a Cu^"^ ion. 

5. Summary 

We have presented the results of Cu NQR and NMR in a 
bicomponent magnet Cu2Fe2Ge40i3 including quantum spin 
dimers and classical spin chains. The results of Cu NQR in 
a reference compound Cu2Sc2Ge40i3 containing only spin 
dimers have also been reported. 

In the paramagnetic state of Cu2Fe2Ge40i3, the *'^Cu nu- 
clear spin-lattice relaxation rate I/Ti nqr has two indepen- 
dent contributions. One is from Fe^^ spins which is indepen- 
dent of temperature, and the other is due to Cu^^ spins ex- 
hibiting a thermally-activated behavior except near T^. The 
temperature dependence of the latter is described by a spin 
gap identical to that in Cu2Sc2Ge40i3. Such a characteristic 
of 1 /Ti NQR suggests that the Fe chains contribute almost in- 
dependently to the magnetic excitation spectrum without af- 
fecting the gapped excitations of the Cu dimers. This is con- 
sistent with the results of the neutron scattering experiments 
and should be ascribed to different energy scales between the 
Cu dimers and the Fe chains. 

There appears an internal field at the Cu site in the anti- 
ferromagnetic state accompanied by ordering of Cu magnetic 
moments. The internal field is nearly parallel to the Z prin- 
cipal axis of the EFG tensor, but is directed away from the 
ordered moment of the Cu atom. For this property of the in- 
ternal field, we introduced a simplified model of the internal 
field based on the fourfold planar coordination of Cu. The 
model revealed that the relation between the directions of the 
internal field, ordered moment and the Z principal axis of the 
EFG tensor is well described by an anisotropic hyperfine in- 
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compatible with a modestly large value of the internal field if 
we take account of a large canting angle of the ordered mo- 
ment from the Z axis as well as a reduction of the moment 
from a fully polarized value. The d(x^ - y^) orbital also pre- 
dicts a value of I/^i^nqr in the high- temperature limit which 
agrees well with that evaluated from the experimental data. 
We therefore conclude that the ground state orbital of a Cu^"^ 
ion in Cu2Fe2Ge40i3 has predominantly a dix^-y^) character. 
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